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ABSTRACT 

Recently, scale-scale correlations have been detected in the distributions of 
quasar's Lja absorption lines and the maps of cosmic temperature fluctuations. 
In this paper, we investigate the scale-scale correlations in galaxy distributions. 
Using samples of mass field given by N-body simulation, we first show that 
the scale-scale correlation of 2-D and 3-D mass distributions at present day is 
capable of breaking the degeneracy between the SCDM (standard cold dark 
matter model) and OCDM (open CDM model) or LCDM (flat CDM model), and 
even show the difference between the OCDM and LCDM. Using biased galaxy 
samples produced in an appropriate bias model, we show that the scale-scale 
correlation of galaxy distribution at zero redshift is still a powerful tool to break 
the degeneracy in the parameter space, including both cosmological and biasing 
parameters. We analyze the scale-scale correlations of the APM bright galaxy 
catalog, and compare it with the mock catalog in the SCDM, OCDM and 
LCDM models. We find that all the spectra of local and non-local scale-scale 
correlations predicted by the OCDM model are in excellent agreement with 
those of the APM-BGC sample, while the SCDM and LCDM mock catalog 
appear to have somewhat weaker scale-scale correlations than the observation. 

Subject headings: cosmology: theory - galaxies: statistics - large scale structure 
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1. Introduction 

A problem in current study of the large scale structures of the universe is the so 
called degeneracy among viable models. That is, structure formation models with different 
cosmological parameters are able to reproduce the same features as the observed large-scale 
structure. For instance, in the cold dark matter (CDM) cosmogony, the standard CDM 
model (SCDM) predicts about the same abundance and two-point correlation functions 
of the present (z = 0) galaxy clusters as the open CDM model (OCDM) and the low 
density flat model (LCDM). This is sometimes called the degeneracy in parameter space 
of Q and cr 8 , where Q is the present value of the cosmological density parameter, and <t 8 
the normalization of the power spectrum on a 8 hr l Mpc scale. Considering the redshift 
evolution of clusters, the degeneracy between the SCDM and OCDM (or LCDM) is broken 
(e.g. Jing & Fang 1994, Bahcall, Fan & Cen 1997.) However, the models of OCDM and 
LCDM remain in degeneracy for the second and lower order statistical features of clusters 
till redshift s as high as z ~ 0.5. 

In terms of the clustering of galaxies, model-degeneracy is more common. For instance, 
the mock catalog of SDSS and 2dF galaxy redshift surveys were generated by a bias 
sampling of the simulated mass fields. The parameters used for the bias sampling were 
selected by fitting the two-point correlation function of the bias sampled galaxies with 
observations (Cole et.al, 1998.) Thus, a majority of the dark matter models are found 
to be able to produce mock catalog which is basically in agreement with the observed 
power spectrum and two-point correlation function of galaxies if the biasing parameters 
are appropriate selected. Namely, to a first approximation, all the models considered are 
degenerate. This is a degenerate in the space of both dark matter and biasing parameters. 

It is largely believed that the non-linear behaviors of cosmic clustering would be 
useful for the degeneracy breaking. In this paper, we will investigate the model-degeneracy 



-4- 



breaking by a non-linear measure, the local scale-scale correlations. The hierarchical 
scenario of structure formation is characterized by a rule that determine how are the small 
scale massive halos related to their parent halos on large scale. The relation between 
parent and daughter halos leads to a local scale-scale correlation. Different models are 
characterized by different relations between massive halos on different scales. Therefore, 
the local scale-scale correlation is sensitive to models. This motivated us to study the 
possibility of breaking the model-degeneracy by the local scale-scale correlations. 

Moreover, the local scale-scale correlations have been detected in 1-D samples like Lj-a 
forests (Pando et al. f998), and 2-D samples from COBE-DMR cosmic temperature map 
on angular scales of fO - 20 degrees (Pando, Vills-Baude & Fang 1998). Namely, the current 
data of large scale structures are ready for detecting scale-scale correlations. In this paper, 
we show that the scale-scale correlations of the APM-BGC galaxies can provide meaningful 
information for degeneracy breaking of cosmological models. 

The outline of this paper is as follows. Based on the discrete wavelet transform (DWT), 
§2 gives a brief description of the basic idea and algorithm for the scale-scale correlations. 
In §3 the scale-scale correlations are computed in simulated samples for three typical CDM 
models. The scale-scale correlation of underlying dark matter and biased galaxies in real 
space as well as redshift space are investigated in detail. The comparison of the scale-scale 
correlations between the APM-BGC galaxies and the model predictions are presented in §4. 
Finally in §5, we briefly summarize the conclusions drawn from our calculations. 



2. Basic formulae of scale-scale correlations 



In the hierarchical clustering scenario, the evolution of massive halos is prescribed 
by an evolutionary tree of halos, such as smaller halos merging into a bigger one. This 
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process leads to correlations between density perturbations on different scales, but 
spatially localized. The correlation can easily be described in phase space, i.e., the 
position-wavevector space. For phase space analysis, a density field p(r) is decomposed 
into perturbations at phase space points (k, x), where k and x denote, respectively, the 
wavevector and position of a volume element d 3 xd 3 k ~ 1 in phase space. The localized 
scale-scale correlation characterizes the correlations between the density perturbations at 
two phase space points with the same spatial coordinate x, but different scale coordinates 
ki and k 2 . 

One way to realize the phase space decomposition is by the discrete wavelet transform 
(DWT) (Fang & Thews 1998). With the DWT decomposition, a 3-D density field p(x), 
x = (xi, x 2 , x 3 ), in the range < Xi < Lj, {i — 1, 2, 3} can be expanded as 

oo 2-j-l 

PW=P + EEW X ) (!) 
j i=o 

V>j,i(x) are the bases of the DWT, j = (ji,j 2 ,h), 1 = (h,h,h) and j t = 0,1,.., 

U = 0, 1, ...2 Ji — 1. Because these bases are orthogonal and complete, the wavelet function 

coefficient (WFC), e^, is computed by 

= J p(x)Vj,i(x)dx. (2) 

From the construction of wavelet bases, we know that indices j and 1 corresponds, 
respectively, to the scale (Li/2 ji ,i = 1,2,3) and position x = {Lili/2 ji ,i = 1,2,3}. 
Sometimes we use wavevector k, which defined as k/27r = {2^/Li,i = 1,2,3}. Therefore, 
ij i represents the amplitude of the density fluctuation localized at point (k, x) of phase 
space. 

In general, density fields in cosmology are assumed to be ergodic, i.e., the average over 
an ensemble is equivalent to the spatial average in one realization. Thus, the ensemble 
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average of perturbations on scale j is 

(q> = X>j,i = °- ( 3 ) 
j 

where the normalized summation notation = rit=i S^o* was used. The last step of 
eq.(3) is due to the fact that the bases of wavelet are admissible, or / ^(x)dx = 0. Thus, 
the variance or the power spectrum with respect to the DWT decomposition is determined 
by (Pando & Fang, 1998) 

Pdwt ® = ^r^(^) = Dw- ( 4 ) 

j 

The local correlation between perturbations on scales of j = (ji, J2, J3) and 
j + 1 = (ji + 1, j 2 + 1, jz + 1) can be measured by 

j+1 

where the [ ]'s denote the integer part of the quantity. Because Lil/2 j = Li2//2- ?+1 , the 
position li at the scale ji is the same as the positions 21 1 and 21 1 + 1 at the scale j\ + 1. 
Therefore, eq.(5) measures the correlation between scales j and j + 1 at the same physical 
point x, i.e. it is localized scale-scale correlation. If perturbations on scales j and j + 1 
are independent, we have (ejCj+i) oc (ej)(ej +1 ) = 0. It is interested to point out that the 
statistic of eq.(5) is of second order. 

The statistical properties of ej ; i with respect to 1 generally are independent from the 
statistics on j. For instance one can have distributions for which ej ; i is Gaussian in its 
one-point distribution with respect to 1 while scale-scale correlated in terms of j (Greiner, 
Lipa & Carruthers, 1995). A simple example is as follows. Considering a distribution of ej^i 
which is Gaussian in 1 at a given scale ji. Namely, higher order correlations of ej^i are zero 
at this scale. Let's consider the perturbation ej 2i i at scale j2- If at a given physical point x, 
ej 2 i is always proportional to ej l5 i at the same point, the distribution of ej 2i i is also Gaussian 
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in the 1 distribution. However, the perturbations on scales ji and j 2 are strongly correlated. 
Obviously such localized scale-scale correlations can only be detected by a scale-space 
decomposition. 

Instead of eq.(5), we will employ a normalized scale-scale correlation defined as 

c p,p = _^j+i e j-[i/2] e j-+i,i 

J V"^ ~p V"^ ~p ' ^ ' 

Z^i e j,[l/2]Z^i +1 e j+l,l 



where p is an integer. In the case of odd p, the denominator would be zero [eq.(3)], and 
therefore, we will use even p only. The statistics Cj^ measures the p-order correlations 
between local density perturbations with wavevectors k/2n = (2 J1 /Li, 2 J2 /L 2 , 2 J3 /L 3 ) and 
k'/2n = (2 jl+1 /L 1 ,2 j2+1 /L 2: 2 js+1 /L 3 ). In the case of p = 2, is actually the correlation 
between the DWT spectra [eq.(4)] on scales j and j + 1, i.e. it is a local spectrum-spectrum 
correlations. 

Generally, we can measure the p-order correlation between perturbations at two phase 
space points (j, 1) and (j + 1, 1 + Al) by 



_ j+i e j,tV2] e j+l,l+Al 
Z^j e j,[l/2]Z^ j+ i e j+l,l 



For p = 2 and Al ^ 0, C?+ 1A1 describes non-local spectrum-spectrum correlations. 



Obviously, C? ,p = 1 and Cj^ A1 = 1 for Gaussian fields. The inequality Cj 3 ' p > 1 and 
Cj+i,Ai > 1 correspond to a positive scale-scale correlation, and Cf' p < 1 and Cj^ p liA1 < 1 
to the negative case. Despite the statistics of eqs. (6) and (7) being higher order for p > 2, 
they are essentially the two-point statistics in phase space. 

In this paper, we will consider only diagonal components of scale-scale correlation, i.e., 
L 1 = L 2 = L 3 , ji = j 2 = j 3 = j and A/ x = Al 2 = A/ 3 = Al. 



- 8- 



3. Scale-scale correlation of simulation samples 

3.1. Samples of N-body simulation 

To demonstrate the scale-scale correlations as a tool of dengeneracy breaking, we first 
calculate the scale-scale correlation of mass fields given by N-body simulation. We used the 
modified version of AP3M code (Couchman, 1991). The cosmological models considered 
here are standard cold dark matter model (SCDM), low density open model (OCDM) and 
low density flat model (LCDM). The model parameters (f2 , T, a 8 ) are taken to be 
(1.0,0.0,0.5,0.55), (0.3,0.0,0.25,0.85) and (0.3,0.7,0.25,0.85) for the SCDM, OCDM and 
LCDM, respectively. 

The simulations were performed by evolving 128 3 particles in a periodic cubical box of 
side 256/i~ 1 Mpc. The soft parameter is taken to be 15% of grid spacing, i.e. 300/2T 1 kpc, 
which keep unchanged in the comoving coordinates, and the force resolution approximates 
100 h~ x kpc for a Plummer force law correspondingly. The initial perturbations were 
generated by the Zeldovich approximation applying to a particle distribution with a 'glass' 
configurations (White 1994). The time steps is 400 for SCDM evolving from z — 15 to the 
present, and 600 for LCDM and OCDM from z = 20. 

These models are degenerate at low redshift z ~ 0, in the sense that they produce 
almost the same abundance and power spectrum of clusters at z — 0. The degeneracies of 
SCDM - OCDM and SCDM - LCDM may be broken if the abundance of clusters at high 
redshifts is considered. However, both LCDM and OCDM predict almost the same redshift 
evolution of cluster abundance until redshift z ~ 0.5. 
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3.2. Scale-scale correlation of 3-D mass distribution 

To probe the scale-scale correlation, we placed randomly 1,000 cubic box of 64 3 h~ 3 
Mpc 3 in the 256 3 /i~ 3 Mpc 3 simulation space. These boxes form an ensemble of our statistics. 
Each 64 3 box is divided into 128 3 grids, and the 128 3 3-D matrix of density distribution is 
then obtained by the triangular shaped cloud (TSC) mass assignment scheme. We subject 
the 3-D matrix of density field by a Daubechies 4 transform, and find the WFCs e^i, where 
j =1,2. ..6, corresponds to scales (64/2 J ) hr 1 Mpc. 

The scale-scale correlation defined by eq.(6) is calculated for the ensemble of 64 3 hr' i 
Mpc 3 boxes in all the three cosmological models. Figure 1 plots the histogram of the 

2 2 

distribution of logCy for j = 2 to 5 in the OCDM model. It shows that the scale-scale 
correlation increases with the j, i.e. the scale-scale correlation on smaller scale is greater 

2 2 1 

than that on larger scales. The distribution of Cj' at j = 2, or ~ 16 h~ Mpc, appears to 

2 2 

close a log- normal distribution. However, the mean value of C 2 ' already departs from one. 
At j = 3, i.e, at the spatial scale ~ 8 h^ 1 Mpc, the confidence level of C 2 ' 2 > 1 is > 95%. 

2 2 

Down to the minimum sturcture resolved here, the distribution of C 5 ' is peaked as high as 
about 30, which indicates that the non-linear evolution of the density field on small scales 
leads to a significant positive scale-scale correlation. 

We now examine the model-dependence of the scale-scale correlation. Since the 
distributions of log C 2 ' 2 for j > 3 are highly non-Gaussian, the mean value and variance of 

2 2 

logCy is not enough to describe their statistical behavior. We introduce the cumulative 
distribution of logC 2 ' 2 , defined by 

F(<1 ° gC < >- epilog CfMlogCf ' (8) 

Figure 2 compares the cumulative distribution F(< logC 2,2 ) of simulation samples in the 
SCDM, OCDM and LCDM models at z = 0. This figure shows clearly that the curve 
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F{< logC 2 ' 2 ) of SCDM is significantly different from those of LCDM and OCDM. That is, 
low redshift (z ~ 0) data are already capable of breaking the degeneracy between SCDM 
and LCDM as well as OCDM if the scale-scale correlation of mass distribution could be 
detected. Moreover, Figure 2 also shows that the curves F(< logCr) for LCDM and 
OCDM are distinguishable from each other at low redshift, despite the difference between 
LCDM and OCDM is small relatively. This results, at least, demonstrates the possibility of 
LCDM-OCDM degeneracy breaking by the DWT scale-scale correlations. 

The power of scale-scale correlation increases from SCDM to LCDM and then OCDM. 
This is consistent with the sequence of non-linear evolution. The SCDM clustering is less 
non-linear than LCDM and OCDM, and the OCDM model shows somewhat more compact 
structures than LCDM (Jing et al. 1995). 

3.3. Scale-scale correlation of 2-D mass distributions 

To measure the scale-scale correlation in 2-D distribution, we make an ensemble of 
samples by randomly placing 1,000 128 x 128 x 16 h~ 3 Mpc 3 slabs in the simulation space. 
A mass distribution in 2-D space of 128 x 128 (h^ 1 Mpc) 2 can be obtained by summing up 
all masses along the dimension of 16 hr 1 Mpc. The density distribution is tabulated by 
TSC mass assignment scheme on 1024 2 grids in the 2-D. Accordingly, the j runs 1,2,. ..,9, 
which corresponds to scales (128/2' ? ) hr 1 Mpc. 

Figure 3 presents the histogram of the log C 2 ' 2 distribution for the 2-D samples in 
the OCDM model. Obviously, the figure shows the similar behavior as 3-D scale-scale 

2 2 

correlation. The power of the 2-D scale-scale correlations Cj' is almost the same as that 
for the 3-D samples but with the slightly smaller values than in the 3-D case. Generally, the 
ordinary two-point correlation strength in 2-D is significantly smaller than in 3-D due to the 
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projection effect. However, the Cy measures an ratio between localized correlations, the 
effect of projection is weaker. Therefore, detecting scale-scale correlations in 2-D samples 
would be as effective as in 3-D cases. 

The model-dependence of Cy 2 in the 2-D samples are shown in Figure 4, which 
compares the cumulative distribution F(< \ogC 2,2 ) between the SCDM, OCDM and 
LCDM models at z = 0. Similar to the 3-D case, the curve of F(< fogC 2 ' 2 ) for the SCDM 
is substantially different from that of the LCDM and OCDM models. The curves for the 
LCDM and the OCDM are also distinguishable. This result does not very sensitively 
depend upon the thickness of the slab. Using the slab of 128 x 128 x 32 h~ 3 Mpc 3 , we 
found very similar results. Thus, to break the degeneracy between the SCDM, LCDM and 
OCDM models, the 2-D data at low redshift (z ~ 0) are as powerful as 3-D samples. 



3.4. Scale-scale correlation of biased galaxy distributions 

Now we turn to the scale-scale correlation of simulation samples of galaxies. In order 
to produce galaxy samples comparable with the mock samples of SDSS and 2dF galaxy 
redshift surveys in Cole et al. (1998), we employ the so-called two-parameter Lagrangian 
bias model to identify the galaxies from N-body simulation mass distributions. Namely, 
the density field is first smoothed by a Gaussian window function of width 3/i -1 Mpc. 
The biased galaxies distribution is then drawn from the smoothed density field p(r) by a 
selection probability defined as 

exv)(au + 3u^ 2 ) if v > 0, 
P(i/)oJ V ' (9) 

exp(az/) if v < 

where the dimensionless variable v is defined by v(v) = 8(r)/a; 5(r) = (p(r) — p)/p the 
density contrast of the density field p(r), and a 2 = (\5g\) its variance. The two parameters a 
and (3 are determined by the least square-fitting of the observed variance of galaxies count in 
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cubic cell of 5 hr l Mpc and 20 hr l Mpc. The adopted values <7 5 = 2.0 and (x 2 o = 0.67 were 
obtained for APM galaxies samples by Baugh and Efstathiou (1994). This two parameter 
bias models are more flexible in matching both of the amplitude and the slope of the galaxy 
two-point correlation function on scales 1 — 10 hr 1 Mpc, though their underlying matter 
distribution may be quite different. That is, all models used to producing the biased galaxy 
samples are degenerate in terms of the low order statistics such as two-point correlation or 
power spectrum of galaxies. 

Proceeding in same way as used in previous subsection, we placed randomly 1,000 cubic 
boxes of 64 3 h~ 3 Mpc 3 in the simulation space. The biased galaxy samples from these 64 3 
boxes form an ensemble for our statistics of the scale-scale correlations. The distributions 

2 2 2 2 

of Cy- for the OCDM model are presented in Figure 5. The Cy shows the similar behavior 
of scale-dependence as the mass distribution, i.e. the scale-scale correlation increases with 
the j. However, the values of C|' 2 are less than their parent mass distributions. 

We also produced biased galaxy samples in redshift space. For simplicity, the redshift 
distortion is applied along one axis of 3D Cartesian coordinates, which is in fact equivalent 
to the assumption of an infinitely distant observer. The 1000 randomly placed 64 3 h~ 3 Mpc 3 
boxes form an ensemble for calculating scale-scale correlations. The results are also shown 
in Figure 5. The scale-scale correlations in redshift space are generally less than those in 
real space, i.e, redshift distortion suppresses the scale-scale correlation. The difference of 
the scale-scale correlation between the redshift and real spaces disappears gradually with 
increasing scales. 

This result is just what we expect. The local scale-scale correlations are sensitive 
to the accuracy of distance or position measurement. A radial distance uncertainty in 
redshift space is caused by the peculiar velocity. Therefore, the scale-scale correlation will 
be significantly contaminated on scales comparable with the uncertainty of radial distance, 
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while the correlation on large scale will not be affected by the redshift distortion. 

To reduce the effect of redshift distortion, one can measure the scale-scale correlation 
in 2-D samples, as 2-D sample is not affected by redshift distortion. As shown in §3.3, 
the power of the scale-scale correlation of 2-D mass distribution is about the same as 3-D 
sample in real space. Therefore, 2-D sample probably is more suitable than 3-D samples to 
detect scale-scale correlations on scales where the redshift distortion is significant. 

3.5. Non-linearity of bias 

An advantage of the statistics Cj' p is free of linear-bias parameters. Linear bias 
parameter b is defined by 

<J(x) ff = W(x) m , (10) 

where 5(x) g and 5{x) m are the density contrasts of galaxies and density dark matter, 
respectively. Since the DWT is linear, one has 

ej,i g = be Jt i m - (11) 

Therefore, we have 

Cf P g = Cf* m (12) 

If the linear bias is introduced by the variance of perturbations, i.e. 

9 = ^j,l rni (13) 

we still have eq.(12). 

Accordingly, one can detect the non-linearity of a bias model by comparing the C?' 2 of 
the biased galaxy samples with the underlying unbiased matter distribution. As have been 

2 2 

seen from Figures 1 and 5, the values of Cj' of the biased galaxies are significantly less 
than their parent mass distributions. Therefore, the bias model eq.(9) is highly non-linear. 
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2 2 

Figure 6 gives the spectra of for both biased and unbiased galaxy samples in the 
OCDM model. Obviously, the unbiased sample has much stronger scale-scale correlation 
than biased galaxies created by eq.(9). 

4. Scale-scale correlation of APM-BGC galaxies 

4.1. Samples of APM-BGC galaxies 

Now, we attempt to measure the scale-scale correlation in the observational galaxy 
samples, and then make comparison with the model predictions. We analyzed the 2-D 
samples of galaxies listed in the APM bright galaxies catalog (Loveday et al. 1996), which 
gives positions, magnitudes and morphological types of 14,681 galaxies brighter than 16 m .44 
over a 4,180 deg° area in 180 Schmidt survey fields of south sky. The completeness is about 
96.3% with a standard deviation 1.9% inferred from carefully checked 12 fields. Therefore, 
it is large and uniform enough for a 2-D DWT analysis. 

In order to carry out a 2-D DWT analysis, we chose two fields SI and S2 from the entire 
survey area as shown in Figure 1 of Fang, Deng & Xia (1998). SI and S2 do not overlap 
from each other. These two fields are selected to be squares on an equal area projection 
of the sky. The equal area projection keeps the surface number density of galaxies on the 
plane being the same as that on the sky. The sizes of SI and S2 have been taken to be as 
large as possible to cover the whole area of the survey. Each field has angular size of about 
37° x 37°. SI and S2 contain 1,095 and 1,055 elliptical and lenticular (EL) galaxies, and 
2,186 and 2,092 spiral (SP) galaxies, respectively. The galaxies in the regions SI and S2 are 
completely independent. 

We divide each square region into 2 10 x 2 10 (1024 2 ) cells, and the angular scale labeled 
by j is 2 10 ~ J in unit of the size of a cell, j runs from 1...10. Using the estimation of the 
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mean depth of the sample given by the luminosity function from the Stromlo-APM redshift 
survey (Loveday et al. 1992), the cell size is found to be about 65 h _1 Kpc, which is fine 
enough to detect statistical features on scales larger than 1 h _1 Mpc. 

In doing statistical analysis of the APM-BGC, the effect of the 1,456 holes drilled 
around big bright objects should be taken into account. For instance, random samples are 
generated in the area with the same drilled holes as that in the original survey. This may 
lead to a underestimation of the scale-scale correlations. We also removed a few galaxies 
placed in the holes drilled. 

4.2. Sales-scale correlations of EL and SP galaxies 

2 2 

Figure 7 plots the Cj ' of the EL and SP galaxies in the fields I and II. The confidence 
level is estimated from 1,000 randomized samples with the same number of galaxies as 
real data. Both randomized and bootstrap resampling give about the same confidence 
level (Fang, Deng & Xia 1998). Figure 6 shows that the randomized distributions are 
also substantially scale-scale correlated for j = 5, 6, i.e. C 2 ' 2 , Cg' 2 > 1. This is because 
the random sample on small scales (j = 5, 6) essentially is a Poisson process which is 
non-Gaussian and thus leads to scale-scale correlation. On larger scales (j < 4), the mean 
number of galaxies in a j cell is high enough, and their one-point distribution of random 
samples is Gaussian. 

Despite the shot noise, Figure 7 shows that the distributions of the EL and SP galaxies 
of the APM bright galaxies catalog are highly scale-scale correlated at j = 5, 6. The 
scale-scale correlations are also significant (> 95% confidence level) on scales j = 3 and 4, 
which corresponds to the scale ~ 8 h _1 Mpc. The EL and SP galaxies basically have about 
the same behavior of the j-dependence of C 2 ' 2 . It increases with the increase of j. This is 
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qualitatively consistent with our simulation results (§3). In addition, ELs generally have 
higher power of Cj' 2 than SPs. This is consistent with the segregation bewteen ELs and 
SPs. ELs concentrate in clusters, having higher clustering, while the field galaxy SPs are 
lower. 

4.3. Mock samples of galaxies 

For a quantitative confrontation between observed and model-predicted scale-scale 
correlations, we use the mock samples for the SDSS redshift survey (Cole et al. 1998). 
The radial selection effect is prescribed by the Schechter form of the luminosity function 
with the parameters drawn from the the Stromlo-APM bright galaxy survey (Loveday, 
et.al, 1992). Therefore, it is appropriate to use these mock samples for demonstrating the 
model-degeneracy breaking by the APM-BGC galaxy scale-scale correlations. 

We consider still the three popular cosmological models of SCDM, LCDM and OCDM 
here (referred as E4, L3S and 03S respectively in Cole et al, 1998). To mimic APM-BGC 
sample, we extract 10 subsamples out of the SDSS mock catalogue by imposing magnitude 
limit of Bj < 16.44 for each model of the SCDM, LCDM and OCDM. The size of the 10 
fields is the same as the fields I and II. The 10 fields are chose to cover the projected survey 
area as large as possible but have less overlap from each others. This selection ensures the 
statistcal independence to a maximum extent. 

The surface number density of galaxies for all the mock samples is found to be 
approximately the same as that obtained for APM-BGC data, the difference between the 
number densities of the real data and mock samples is less than 3%. The two-point angular 
correlation functions of the mock samples in the three cosmological models are also found 
to be in good agreement with that of APM-BGC. That is, in terms of the first (number 
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density) and the second order (two-point correlation function) statistics, the mock samples 
of the SCDM, LCDM and OCDM are degenerate. 

4.4. The APM-BGC scale-scale correlation and cosmological models 

Since the mock samples do not contain the information of morphology, we calculated 
the scale-scale correlations of APM-BGC galaxies consisting of both ELs and SPs in the 
fields I and II. The results are shown in Figure 8. The 95% confidence level is estimated 
from 1,000 Poisson-distributed samples with the same number of galaxies as the real data. 
The scale-scale correlations of the 10 mock samples in the SCDM, LCDM and OCDM 
models are displayed by the scatter symbols in Figure 8 repectively. 

2 2 

Figure 8 shows first that the fields I and II give almost the same spectrum of Cj' vs. 
j. The data are qualified for model discrimination. At the first glance, all mock samples 

2 2 

of the three CDM cosmogony models show the same shape of the Cj' spectrum, and are 

2 2 

generally consistent with the observed Cj' . However, SCDM, LCDM and OCDM predicted 

2 2 

different Cj' at j = 5. The OCDM is in good agreement with the entire spectrum of 
APM-BGC's Cf 2 from j = 2 to 5, i.e. scales ~ 2 to 16 h" 1 Mpc. However, the values of 
C;/ of SCDM and LCDM are somewhat lower than the observation at 90% confidence 
level. Considering the possible underestimate of the APM-BGC's scale scale correlation 
due to the holes drilled around big bright objects, the difference between observed results 
and the SCDM and LCDM prediction might be more significant than 90%. Thus, the test 
of scale-scale correlation seems to be in favor of the OCDM model. 

It might be too early to say this conclusion is concrete, as the data is still not perfect 
enough. However, this result showed that the scale-scale correlation is effective and feasible 
as a tool of degeneracy breaking. The correlations are able to reveal not only the difference 
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between the SCDM and OCDM, but also the LCDM and OCDM. 

All Cy shown in Figs. 1-8 are calculated by taking Al = [see eq.(7)]. It describes 
the correlations between j and j + 1 perturbations at the same physical coordinates, i.e. 
it is local scale-scale correlation. We also calculated non-local scale-scale correlations by 
taking Al — 1, 2, and 3. It describes the correlation between perturbations on scale j 
at position I and on scale j + 1 at I ± Al. Figure 9 compares the non-local correlations 
c ?,aj=i,2,3 of APM-BGC data and the OCDM mock sample. It shows that the OCDM result 
is in excellent agreement with observed data. 

5. Discussions and conclusions 

The ordinary two-point correlation function is two-point statistics in physical space. 
The scale-scale correlations essentially are two-point statistics both in physical space and 
wavenumber space, i.e., in a position-scale phase space. Moreover, the scale-scale correlation 
describes the feature of the non-linear evolutions of the hierarchical clustering. Therefore, 
the scale-scale correlation is probably the most basic measure of the cosmic clustering in 
the phase space. 

The abundance and two-point correlation functions of galaxies and clusters are unlikely 
to discriminate among the CDM family of models, which are degenerate in cosmological 
parameter space. In this paper, using the mass distribution generated by N-body 
simulation, we showed that the scale-scale correlation of mass distributions is sensitive to 
the cosmological parameters. The scale-scale correlation of biased galaxy distribution is 
sensitive to cosmological and biasing parameters. Therefore, scale-scale correlation of galaxy 
distributions provides us a tool of discriminating among various cosmological models, which 
are degenerate in the parameter space, including both cosmological and bias parameters. 
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Using APM-BGC data, we found that the j spectrum of local and non-local scale-scale 
correlation can generally be matched by the mock galaxy samples of the CDM cosmogonies. 
Amongst the models, the best fitting is given by the OCDM model. Both of SCDM model 
and LCDM model show somewhat weaker scale-scale correlations than OCDM on small 
scales. We can conclude that the OCDM plus the biasing algorithms (9) is favored by the 
current detection of APM-BGC's scale-scale correlations. This showed that the scale-scale 
correlation is a capable tool of degeneracy breaking. 

We found that the scale-scale correlation in 2-D samples is comparable with that in 
3-D. Because the redshift distortion contaminates the scale-scale correlation in 3-D, and 
meanwhile the 2-D samples is free of redshift distortion, the scale-scale correlation measured 
in 3-D redshift galaxy survey will be lower than their 2-D scale-scale correlation. Therefore, 
the difference between 2-D and 3-D may provide us a new way to study the large scale 
velocity field. 

We thank Dr. Wolung Lee for useful discussion. LLF acknowledge support from the 
National Science Foundation of China(NSFC) and World Laboratory scholarship. ZGD is 
supported by NSFC and exchange program between NSFC and DFG. This project was 
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2 2 

Fig. 1. — The histogram of the scale-scale correlations, Cj' , of mass distribution for OCDM 
model. C 2 ' 2 is computed for 1000 randomly throwing box of 64 hr 1 Mpc on a side. 

2 2 

Fig. 2. — The cumulative distribution of logC 5 ' of mass distributions for models of SCDM, 
OCDM and LCDM. 

Fig. 3. — The histogram of the scale-scale correlations, C 2 ' , of 2-D mass distribution for 
OCDM model. C 2 ' 2 is computed for 1000 randomly throwing box of 64 h^ 1 Mpc on a side. 

m 2 2 

Fig. 4. — The cumulative distribution of logC/ of 2-D mass distributions for models of 
SCDM, OCDM and LCDM. For the same type of line, the left and right corresponds to the 
slabs of thinkness 32 h' 1 Mpc and 16 h~ l Mpc, respectively. 

Fig. 5. — The histogram distribution of the scale-scale correlations, C 2 ' 2 , of OCDM biased 
galaxy distributions in both real space and redshift space (filled histogram). 

Fig. 6. — The spectrum of 3-D scale-scale correlation in the range with 90% confidence 
level both for the underlying matter distribution (dark gray band) and biased galaxies (light 
gray band) in OCDM model. The solid lines are given by ensemble average of log C 2 ' 2 . In 
addition, the scale-scale correlation for Poisson-distributed galaxy samples(black band) is 
also plotted. 

Fig. 7. — The 2-D local (A/ = 0) scale-scale correlation of EL and SP samples in areas I 
and II of APM-BGC survey. Gray band shows the range with 95% confidence level for 1,000 
Poisson-distributed galaxy samples. 

Fig. 8. — The 2-D local (A/ = 0) scale-scale correlation of ELs plus SPs in areas I and II 
of APM-BGC galaxy sample (solid lines). The square represents the value obtained from 
one of the 10 mock samples of the indicated model (SCDM, OCDM and LCDM) and are 
displayed in the range with 90% confidence level. The gray bands show the range with 95% 
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confidence level for Poisson-distributed galaxy samples. 

Fig. 9. — The 2-D non-local scale-scale correlation with A/ = 1, 2, and 3 in OCDM. The 
solid line, squares and gray band are the same meaning as those in Fig. 8. 
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